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HIGHLIGHTS 


Model developed to describe the energy balance during willow torrefaction. 

First stage is exothermic releasing 40-280 kj/kg. 

► Second stage is either exothermic or endothermic. 

► Higher torrefaction temperatures result in increased reaction rate and heat release rate. 
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Based on the evolution of volatile and solid products predicted by a previous model for willow torrefac¬ 
tion (Bates and Ghoniem, 2012) a thermochemical model has been developed to describe their thermal, 
chemical, and physical properties as well as the rates of heat release. The first stage of torrefaction, asso¬ 
ciated with hemicellulose decomposition, is exothermic releasing between 40 and 280 kj/kg init iai- The sec¬ 
ond stage is associated with the decomposition of the remaining lignocellulosic components, completes 
over a longer period, and is predicted to be either endothermic or exothermic depending on the temper¬ 
ature and assumed solid properties. Cumulative heat release increases with the degree of torrefaction 
quantified by the mass loss. The rate of mass loss and rate of heat release increase with higher temper¬ 
atures. The higher heating value of volatiles produced during torrefaction was estimated to be between 
4.4 and 16 MJ/kg increasing with the level of mass loss. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Raw biomass faces utilization challenges due to its low bulk en¬ 
ergy density, hydrophobicity, propensity to decay during storage, 
and high grinding energy requirement. Torrefaction has been pro¬ 
posed as thermal pretreatment to improve these characteristics 
(Arias et al„ 2008; Bergman et al., 2005; Ohliger et al., 2012). Tor- 
refaction is a mild pyrolysis occurring between 200 and 300 °C for 
a residence time between several minutes to about an hour which 
results in partial devolatilization (0-60 wt.%) of the solid biomass 
(Prins, 2005). Pyrolysis broadly refers to partial to complete devol¬ 
atilization (70-90 wt.%) occurring over a wide temperature range 
(200-1000 °C) though typically between 400 and 700 °C (Neves 
et al., 2011). 

Multi-scale (e.g. particle, reactor, flowsheet) models of torrefac¬ 
tion require quantification of the physical and chemical properties 
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of torrefaction products in order to describe the process energy 
balance. At the particle scale, the heat of torrefaction may act as 
a local source or sink of thermal energy and therefore affect tem¬ 
perature and conversion profiles significantly. As Turner et al. 
(2010) note, the heat release is a crucial mechanism through which 
chemical reaction and heat transfer phenomena are coupled at the 
particle scale. Another important level of coupling occurs between 
the particle temperature and the solid physical properties, such as 
specific heat, density, and conductivity, which determine the rate 
at which heat diffuses through the particle. Description of the evo¬ 
lution of these properties requires accurate kinetics and thermo¬ 
chemical models. 

At the reactor scale, an accurate estimate of the reaction enthal¬ 
py is necessary for process design and control. For example, pro¬ 
cess conditions which lead to excessive heat release could 
potentially cause runaway thermal reactions in the absence of ac¬ 
tive thermal management. 

For the development of flowsheet models, it is necessary to 
quantify the thermal energy input required during torrefaction to 
assess the overall process efficiency. In particular, the feasibility 
of autothermal torrefaction as originally described by (Bergman 
et al., 2005), which combusts the volatile products released during 
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torrefaction to provide heat for the process itself, depends on the 
accurate estimation of their composition and heating value. 

The present work describes a model for the energy balance and 
thermochemistry of torrefaction based on a previous kinetics mod¬ 
el (Bates and Ghoniem, 2012) for the evolution of volatile matter 
during willow torrefaction. The parameters of the two models 
are linked. While the volatile composition parameters from are 
specific to willow torrefaction, this work outlines a general frame¬ 
work for relating these feedstock-specific parameters to the energy 
balance of the reactions. 

A variety of experimental studies have examined the effect of 
torrefaction conditions on the mass and energy yield of the solid 
product (Almeida et al., 2010; Arias et al„ 2008; Bridgeman et al., 
2008; Medic et al., 2012). Lee et al. (2012) and Medic et al. 
(2012) performed multiple regressions relating process parameters 
(such as initial moisture content, temperature, and residence time) 
to the solid product energy and mass yield. (Almeida et al., 2010) 
showed that the solid mass loss could be used as a quantitative 
indicator for the extent of torrefaction. Several product character¬ 
istics including solid energy yield, heating value, and chemical 
composition were linearly regressed against the solid mass loss. 
The validity of these empirical regressions beyond the thermal 
conditions tested is unknown. Additionally they cannot predict 
the energy content of liberated volatiles or the reaction thermo¬ 
chemistry. The latter are the objectives of this work. 


1.1. Previous estimates for the heat of torrefaction 

Estimates for the heat of torrefaction and pyrolysis as measured 
or deduced by a variety methods are summarized in Table 1. The 
main methods include (1) estimation from the heating value of 
the products, (2) differential scanning calorimetry (DSC), (3) fitting 
of coupled single particle models with experimental data and, (4) 
measurement of heat release using custom reactor systems. 

Experimental attempts to assess the heat of willow and beech 
torrefaction through measuring (or estimating) the heating values 
of the initial biomass, solid, and volatile products have been unable 
to precisely classify it as either exothermic or endothermic (Prins, 
2005; van der Stelt, 2011). Overall heat of torrefaction estimates 
based on this method vary widely and from 87-124 (±449) kj/kg 
willow (Prins, 2005) and -1200 to 1500 kj/kg beech (van der Stelt, 
2011). Significant uncertainties arise from the measurement/esti¬ 
mation of the heating value of solid products and assumptions 
regarding leaked products. 


van der Stelt (2011) also used a fixed bed reactor setup to tor¬ 
refy a single 10 cm long, 28 mm diameter cylindrical beech particle 
in order to deduce the enthalpy of reaction. Six thermocouples en¬ 
abled measurement of the particle temperature at positions be¬ 
tween the centerline and surface. Based on the intraparticle 
temperature profile and a simplified single particle model, the 
reaction heat was estimated to be between 0 and -220 kj/kg beech 
for torrefaction temperatures between 200 and 280 °C. 

Ohliger et al. (2012) utilized a continuous screw reactor system 
to measure heat release during the torrefaction of beech wood. 
Though sensitive to the assumed physical properties of the bio¬ 
mass such as the specific heat capacity, it was concluded that the 
heat of torrefaction was between -199 and 148 kj/kg beech. 

1.2. Previous estimates for the heat of pyrolysis 

Pyrolysis refers to devolatilization occurring over a broader 
range of thermal conditions than torrefaction; however, given the 
greater number of investigations examining the heat of wood pyro¬ 
lysis, it is also worth summarizing these efforts. Literature estimates 
for the heat of wood pyrolysis range widely from -2300 kj/kg 
biomass (exothermic) to 418 kj/kg biomass (endothermic) (Turner 
et al„ 2010). The differences have been attributed to the small 
sample sizes (ranging from milligrams to grams), the presence of 
impurities, and experimental conditions (Turner et al., 2010). 

During differential scanning calorimetry (DSC) wood pyrolysis 
experiments ramping from 100 up to 500 °C at 10 °C min 1 (Rath 
et al„ 2003) found wide variations in the heat of primary wood 
pyrolysis depending on the sample mass, feedstock, and condi¬ 
tions. Observing an apparent shift from endothermic to exothermic 
behavior, it was concluded that the overall heat of pyrolysis de¬ 
pended on the competition between an exothermic char formation 
process (-3525 to -3827 kj/kg c h ar ) and an endothermic (936 to 
1277 kJ/kgg as+vo i a tiies) volatile formation process. The overall heat 
of pyrolysis expressed as a linear superposition of an exothermic 
char formation and endothermic volatile formation weighted by 
mass fractions (obtained experimentally or from kinetic expres¬ 
sions) (Rath et al., 2003): 

AH tota , = AH exo y char + AH end0 (l - y char ) (1) 

where AH tota i has units of kj kg b iomass~\ ychar is a product mass frac¬ 
tion with units of (kg char kg biomass _1 ), AH exo and AH endo are the 
apparent exothermic and endothermic heats of pyrolysis (see 
Table 1 ) for the reported values. 


Table 1 

Experimentally measured or deduced enthalpy of reaction for torrefaction and pyrolysis (-) exothermic, (+) endothermic. 


Enthalpy of reaction 

Temperature 

(°cf 

Feedstock 

Method 


Source 

87 ± 449 kj/ltgbiomass 

250 

Willow 

ASTM bomb calorimetry 


Prins (2005) 

150 ± 1350 kj/kgbiomass 

230-280 

Beech 

Estimated through analysis of products and reactant 

van der Stelt (2011) 

0 to -200 kJ/kg biom a SS 

200-280 

Beech 

Deduced from experimental data wit 
model 

h simplified single particle 

van der Stelt (2011) 

-199 to 148 kj/kgbjomass 

270-300 

Beech 

Deduced from experimental data wit 

h continuous screw reactor 

Ohliger et al. (2012) 

936-1277 kj/kggaa+voiatiles 
—3525 to -3827 kj/kg ch ar 

100-500 

Spruce, Beech 

Differential scanning calorimetry 


Rath et al. (2003) 

-293 to -1673 kj/kg mass i„ ss 

275-470 

Beech 

Deduced from experimental data wit 

h single particle model 

Roberts and Clough 
(1963) 

200.8 kj/kgmass loss 

470 

Beech 

Deduced from experimental data wit 

h single particle model 

Rung and Kalelkar 
(1973) 

255 to -20 kj/kg ma ss loss 

300-600 

Wood sawdust 

Deduced from experimental data wit 

h single particle model 

Koufopanos et al. 

(1991) 

25 kj/kgchar,tar^as 

200-850 

(Various) 

Deduced from experimental data wit 

h single particle model 

Haseli etal. (2011) 

-55.3 to -176 kj/kgbiomass 

100-600 

Pine, oak 
sawdust 

Deduced from experimental data wit 

h model of packed sawdust 

Strezov et al. (2004) 
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Various experimental investigations have deduced the heat of 
pyrolysis through fitting of predictions from coupled single particle 
models with observed particle temperature profiles beginning with 
the work of Roberts and Clough (1963). Such methods lead to vary¬ 
ing estimates ranging from endothermic to exothermic (see 
Table 1 ). One issue which contributes to the these discrepancies 
was originally highlighted by Rung and Kalelkar (1973) - and re¬ 
cently corroborated by Haseli et al. (2011) - is related to whether 
the temperature-dependence of the heat of pyrolysis are accounted 
for in the single particle models. Based on the fitting, they con¬ 
cluded that the standard enthalpy of pyrolysis (at a reference 
temperature) is endothermic (25 kj/kg); however, at elevated tem¬ 
peratures, the enthalpy of pyrolysis can be exothermic due to dif¬ 
ferences in the heat capacities of products and reactants (Haseli 
etal., 2011). 

In reviewing these efforts, it is clear that a thermochemical 
model for torrefaction should account for (1) the dependence of 
the heat of torrefaction on the competitive reaction pathways 
(e.g. char, volatile yield), and (2) the temperature dependence of 
the enthalpy of reactions. 

2. Mathematical model 

2.1. Approach 

The present thermochemical model is coupled to an existing 
kinetics model (Bates and Ghoniem, 2012) which describes the 
evolution of solid and volatile products during willow torrefaction 
(summarized in Section 2.2). Based on the detailed composition 
information provided by this kinetics model, the temperature- 
dependent properties of the pseudo components (e.g., enthalpy, 
energy yield) are estimated (Sections 2.3 and 2.4). Then, equations 
describing the energy balance (heat release) in terms of the reac¬ 
tion kinetics are presented (Section 2.5). In summary, the present 
model describes: 

(1) The temperature-dependent thermodynamic properties of 
all pseudocomponents ( A , B, C, VI, V2). 

(2) The evolution of the solid and volatile products, heating val¬ 
ues, and energy yields. 

(3) The rates of (and cumulative) heat release. 

2.2. Solid mass loss kinetics and composition 

The kinetics model developed in (Bates and Ghoniem, 2012) 
used a two-step, first-order mechanism to describe the solid mass 
loss kinetics. The mechanism was originally shown by (Prins, 2005) 
to satisfactorily describe willow mass loss curves during torrefac¬ 
tion. In this mechanism, the torrefaction products are lumped into 
five pseudo-components: 



where the solid phase pseudo components (A, B, C) represent raw 
biomass (A), an intermediate solid (B), and char (C), volatile prod¬ 
ucts are represented by VI and V2, and k,, k 2 , k vi , k^ represent 
Arrhenius kinetic parameters. The solid mass yield Y s and mass loss 
Vioss are obtained through integration of Eq. (1) through Eq. (5) in 
(Bates and Ghoniem, 2012): 
y _ m A -m B - me 

s ~ m 0 U 


Y loss = mV '+ o mV2 = 1 _ Vs (3) 

where m x is the mass of pseudo-component X (A, B, C, VI, V2), and 
m 0 is the initial dry mass of the solid. 

The elemental composition of A is a model input and deter¬ 
mined from the ultimate analysis (dry basis) of the raw biomass 
it represents. For this work, the elemental composition of raw wil¬ 
low from (Prins, 2005) was assumed: Y jA = 0.472, 0.061, 0.451, 
0.003, 0.013 j = C, H, O, N, Ash, respectively. Eq. (25) and (26) in 
(Bates and Ghoniem, 2012) define the composition of B and C 
(Yj.b and Y jA ) in terms of the kinetic rate parameters and composi¬ 
tions of A, VI, and V2. The composition of the average solid prod¬ 
uct, Yjf is calculated from the integration of Eq. (18) (Bates and 
Ghoniem, 2012). Volatile pseudo components VI and V2 are de¬ 
scribed by the fixed compositions shown in Table 2 of (Bates and 
Ghoniem, 2012). 

The characteristics of this two-stage kinetics mechanism are re¬ 
lated to the thermal stability of the components of the biomass. 
Aside from ash, extractives, and moisture, lignocellulosic biomass 
is comprised of three components including hemicellulose, cellu¬ 
lose and lignin. It is suggested that the more rapid first stage of tor- 
refaction can be attributed to the decomposition of the 
hemicellulose while the slower second stage is explained by the 
decomposition of the more thermally stable cellulose fraction. 
(Bates and Ghoniem, 2012; Chen and Kuo, 2011; Prins, 2005; van 
der Stelt, 2011 ). Since ash and lignocellulosic composition varies 
between feedstock type (woody, herbaceous) and species their 
reactivities differ significantly (Bridgeman et al., 2008). The present 
model utilizes kinetics and volatile composition parameters spe¬ 
cific to willow (a hardwood); so cross validation of certain model 
predictions with published data from other lignocellulosic feed¬ 
stocks was carried out. 

2.3. Thermodynamic property estimation 

In order to model the reaction thermochemistry, estimates of 
the total enthalpy for the nine volatile matter species and three so¬ 
lid pseudocomponents were necessary. The total enthalpy at tem¬ 
perature T is given by the summation of the formation enthalpy 
and sensible enthalpy: 

Hi(T) = Hf t + j* c p ,i(T)dT (4) 

where H° f i is the standard heat of formation of species i in J/kg and 
c Pii is the specific heat capacity of component i in J kg -1 K 1 , T is the 
reaction temperature in Kelvin, and T 0 is the standard temperature 
(298.15 K). 

2.3.1. Volatile products (VI, V2) 

VI and V2 were modeled as ideal gas mixtures of nine chemical 
components including acetic acid, water, carbon dioxide, metha¬ 
nol, lactic acid, furfural, hydroxyacetone, carbon dioxide, and car¬ 
bon monoxide. Property data for the heating value, enthalpy of 
formation, and specific heat capacity for these chemical species 
was combined from various sources and summarized in Table Al. 

2.3.2. Solid reactant and products (A, B, C) 

2.3.2.1, Heating value. Several validated correlations have been 
developed to predict the heating values of coal and biomass solid 
fuels based on their elemental composition as defined through ulti¬ 
mate analysis. To the best knowledge of the authors, no correla¬ 
tions have been specifically developed or validated for torrefied 
biomass. The accuracy of the correlations as measured by their 
standard error of prediction (SEP) when compared to a set of bomb 
calorimetry measurements is between 337 and 1023 kj/kg (Friedl 
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Table 2 

HHV correlations based on the ultimate analysis of a solid fuel on a dry basis. 


Number Name of correlation 

Correlation (HHV, J kg- 1 ) 


Source: 


(1) Bolt 


HHV = 1000[351.69(y c ) + 1162.46(y„) - 110.95(y o ) + 

104.67(y s ) - 62.8(y N )] 

Boie (1953) 


(2) Frie 

(3) IGT 

dl et al. 

HHV = 10 6 [35.45(y c ) 2 - 23.15(y c ) - 223.35(y H ) + 512(y c )(y H ) + 13.1(1*,) + 20.5875] 
HHV = 1000[341.7(y c ) + 1322.1(y H ) - 119.8(y 0 + Y N ) - 15.3(Va*)] 

Friedl et al. (2005) 

Institute of Gas Technology (1978) 

Table Al 

Specific heat 0 kg 1 K 
Sources: (Goose et al., 

~’), specific standard enthalpy of formation (J kg ’), and molecular weight for 
2012; National Institute of Standards and Technology, 2011). 

main volatile components (in gaseous phase). VI, 

V2 properties calculated. 

Chemical name 

Formula 

Molecular weight (g mol J ) 

Cp&s (J kg 1 K ’) 

f^akg- 1 ) 

HHV (MJ kg -1 ) 

Acetic acid 

c 2 h 4 o 2 

60.052 

435.3 +2.18T 

-7,197,978 

14.6 

Water 

h 2 o 

18.01528 

1670 + 0.64T 

-13,423,383 

0 

Formic acid 

ch 2 o 2 

46.0257 

474.5+ 1.753T 

-8,225,187 

5.55 

Methanol 

ch 4 0 

32.04186 

601 +2.521T 

-6,271,172 

22.7 

Lactic acid 

c 3 h 6 o 3 

90.07794 

481.6+ 2.514T 

-6,775,000 

15.1 

Furfural 

C5H402 

96.08406 

331.09+ 2.367T 

-1,556,970 

24.4 

Hydroxyacetone 

c 3 h 6 o 2 

74.07854 

494.2 + 2.338T 

-4,833,775 

21.9 

Carbon dioxide 

co 2 

44.0095 

978.4 + 0.2T 

-8,941,479 

0 

Carbon monoxide 

CO 

28.0101 

767 + 0.4T 

-3,946,077 

10.1 


Co. 3 8 8 Hi.9, 

37O1229 26.24 

1162 +0.961 T 

-10,345,151 

4.43 

V2 

Ci 3 62H 3 .9. 

440l. 543 45 

601 + 2.262T 

-7,078,295 

16.3 


et al., 2005) or an average absolute error between 2.6% and 6.7% 
(Sheng and Azevedo, 2005). Based on their performance in litera¬ 
ture evaluations, three equations (shown in Table 2) to calculate 
the HHV in J/kg on a dry basis were selected and compared for 
modeling. Y, is the elemental composition expressed as the mass 
fraction of element j (C, H, 0, S, N, Ash) contained in the dry fuel. 
HHV Eq. (1) in Table 2, the well-known Boie correlation, was devel¬ 
oped originally for coals but is commonly applied to biomass sam¬ 
ples as well. HHV Eq. (2) in Table 2, the Friedl et al. correlation, was 
derived from ordinary least squares regression and partial least 
squares regression fitting of 122 biomass samples. The Friedl 
et al. correlation was fitted/validated for biomass samples with 
higher heating values between 16 and 26 MJ/kg and carbon con¬ 
tent between 42% and 50%. HHV Eq. (3) from Table 2, the IGT cor¬ 
relation, was also originally developed for coal but has since been 
applied and validated for refuse derived fuel (Sheng and Azevedo, 
2005). 

23.2.2. Enthalpy of formation. The enthalpy of formation of solid 
fuels can be estimated based on their heating value. Complete 
combustion of a dry biomass fuel represented by the chemical for¬ 
mula CH x O y can be expressed as: 

CH*O y (s) + ( 2 + ^ 2 ~ y ) 0 2 (g) - *H 2 0(1) + C0 2 (g) (5) 

where the stoichiometric coefficients x and y can be determined 
from the fuel’s elemental composition, 


y H MW c 

(6) 

Y^MW^ 

Y 0 MW C 

(7) 

Y c MW h 


where MW C , MW H are the atomic weights (kg mor 1 ) of carbon and 
hydrogen, respectively. The relationship between the HHV (re¬ 
ported as a positive value) and the enthalpy of formation of reac¬ 
tants and combustion products is: 

-HHV = £ v,% - (8) 

products reactants 

where, HHV is the higher heating value (J mor 1 ), H^is the standard 
enthalpy of formation (J mor 1 ) for the ith product or reactant, and 


Vj is the number of mols product or reactant per mol of solid fuel 
reactant. After substitution and rearrangement, the specific stan¬ 
dard enthalpy of formation (J/kg fuel ) can be defined in terms of its 
higher heating value (J/kgtuei) and elemental composition: 


(9) 


23.2.3. Specific heat capacity. The temperature and composition 
dependent specific heat capacities (J kg -1 K ') of the solid pseudo 
components were modeled according to (Merrick, 1983): 



(10) 


where R is the ideal gas constant 8.314J mor 1 K ’and the average 
molecular weight MW X of the solid pseudo component X (A, B, C ) is 
calculated from its respective elemental composition, 

\/MW x = J2Y DAFj /MWj (11) 


where, Y DAF j, j = 1...,5 represent the mass fractions of carbon 
hydrogen, oxygen, nitrogen, and sulfur on a dry, ash-free basis 

and MWj, j= 1.5 represents the respective atomic weights 

(kg mor 1 ) of the elements. 

2.4. Energy yield 


2.43. Solid product energy yield 

The solid energy yield (in terms of the HHV) of the torrefied so¬ 
lid product is an important torrefaction parameter and various 
published experimental results exist for a variety of feedstocks 
and torrefaction conditions. The definition of the solid energy yield, 
ris is: 

Vs = YsRh hv (12) 

where an energy densification ratio, R H hv is defined: 


(13) 
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Table 3 

Feedstocks, temperature ranges, and methodology for determining torrefaction solid product energy yield from she different publications. 


Feedstock Temperature range (°C) Residence time (min) 


Method of measuring heating 


value 


Prins (2005) Willow 

Almeida et al. (2010) Eucalyptus 

Arias et al. (2008) Eucalyptus 

Bergman and Kiel (2005) Willow 

Medic et al. (2012) Com stover 

Ohliger et al. (2012) Beech 


250-300 17-38 

180-280 60-300 

280 15-180 

280 17.5 

200-300 20-31.4 

270-300 15-40 


ASTM bomb calorimetry 
ASTM bomb calorimetry 
Not specified 
ASTM bomb calorimetry 
Ultimate analysis with correlation 
ASTM bomb calorimetry 


HHV/is the heating value of the torrefied solid product in J kg \ and 
HHV 0 is the initial heating value of the solid feed in J kg ’. 


2.4.2. Volatile product energy yield 

The energy yield of the total volatiles is defined: 


tfvtot — Ylos 


HHVytot 

HHV 0 


(14) 


where HHV vto t is the average heating value of the cumulative vola¬ 
tiles J kg 1 


2.5. Heat release 


The heat release rate for the overall reaction as well as that 
attributed to each stage is calculated from the thermochemical 
properties and reaction rates. The general equation for the rate of 
heat release under isothermal conditions at temperature T can be 
summarized: 

#+£«£*-<> o=) 

where q r is the cumulative heat release in J; H x (X = 1,2,3... 5) is the 
enthalpy of component X at temperature T; m x (X = 1,2,3... 5) rep¬ 
resents the mass of component X. The index X refers to the pseudo¬ 
components A, B, C, VI, V2, respectively. This differential equation 
can be integrated simultaneously with the other mass loss kinetics 
equations to compute the cumulative heat release rate. 

The instantaneous product distribution (at temperature T) is de¬ 
fined by the solid mass loss kinetics mechanism which can be de¬ 
scribed by a two stage first-order decomposition: 

A —> /IB + vVl (16) 

B^yC+ZV 2 (17) 

where fi, v, y, £, the instantaneous fractional yields, were previously 
defined (Bates and Ghoniem, 2012). Therefore an enthalpy of reac¬ 
tion for stages one, two, and at completion under at temperature T 
can be defined: 

AH r> j = /jHb + vHy i — Ha (18) 


AH r ,2 = yH c + cH V2 - H b 


(19) 


AH r ,final = PyHc + vHvi + [1cH V2 - H a (20) 

where, H x is the total enthalpy of pseudocomponent (X = A, B, C, VI, 
1/2) in J kg -1 . The enthalpies of reaction are defined per kilogram of 
reactant; so for step one, AH ril is defined per kg of A, A H r2 is defined 
per kg of B and A H rJina , is defined per kg of biomass. At the comple¬ 
tion of both steps, the products consist entirely of C, VI, and V2. 


3. Results and discussion 

First the energy yield and predictions of the model are com¬ 
pared with experimental results in Sections 3.1 and 3.2 for the so¬ 
lid and volatile product, respectively. Then the model predictions 


for the magnitudes and rates of heat release are discussed in Sec¬ 
tion 3.3. Lastly, a brief evaluation of the model and sources of 
uncertainty are described in Section 3.4. 

3.1. Solid product energy yield and densification 

The model predictions of the solid energy yield ( rj s ) versus mass 
loss (Tioss) ate compared with published experimental data in 
Fig. la. Fig. lb illustrates the energy densification ratio (Rhhv) ver¬ 
sus mass loss (Yi oss ). In Figs, la and b and 2 the model input condi¬ 
tions representative of willow torrefaction (T in i tia i 200 °C, 
Tfinai = 300 °C), and a heating rate of 5.88 °C/min) are assumed. 

The data used to validate the present model was gathered under 
a variety of experimental conditions, feedstocks, and methods for 
measuring the energy yield, which are summarized in Table 3. Sev¬ 
eral of the studies, (Almeida et al„ 2010; Bergman and Kiel, 2005; 
Ohliger et al„ 2012; Prins, 2005), used bomb calorimetry to deter¬ 
mine the energy yield. As highlighted by (Bergman et al., 2005), 
bomb calorimetry has a ±240 kj kg -1 uncertainty in heating value 
determination, therefore the relative percentage uncertainty in 
the higher heating value is ±1.5% and the resulting uncertainty in 
the estimated energy yield is approximately ±3% . 

(Medic et al., 2012) calculated the heating value of the raw and 
torrefied corn stover by measuring the product ultimate analysis 
and applying a HFIV correlation. Elemental analysis has an error 
of ±0.3 wt.% (Friedl et al., 2005). The correlation employed by Med¬ 
ic et al. (2012) accurately predicted the heating value of 90% of 
tested biomass samples within ±5% (Sheng and Azevedo, 2005). 

The data from (Almeida et al., 2010) represents the linear corre¬ 
lation they developed to describe the energy yield based on torr¬ 
efied eucalyptus versus its mass loss: 

JJ S = (1-0.6Y |OSS ) (21) 

The equation suggests that a 10% mass loss results in a 6% loss in 
energy yield. (Arias et al., 2008) did not report the method for deter¬ 
mining the heating value or the uncertainty associated with their 
proposed correlation. 

The model results illustrate a trend observed in all the experi¬ 
mental data: the heating value of the torrefied willow increases 
monotonically as mass loss increases. Applying the Boie and Friedl 
et al. HHV correlations shows good agreement in the predictions of 
solid energy yield. At relatively lower levels of mass loss (0-15%), 
there are very small differences between the predictions. The mod¬ 
el predictions between 10% and 16% mass loss show the solid 
energy yield to be within 1% of the two experimental measure¬ 
ments made by Bergman & Kiel and Prins and within 3% of the 
points from (Medic et al., 2012). As the mass loss increases to 
greater than 25%, the model energy yield agrees with the data 
within the experimental uncertainty (±3%), but appears slightly 
high. For example, after 33% mass loss, the model predicts an en¬ 
ergy densification ratio between 1.25 and 1.24 (for HHV Eqs. (2) 
and (1), respectively) which is slightly above Prins’ experimental 
measurement of 1.19. At this high/intermediate range of mass loss 
(25-35%) It is difficult to conclude whether the model significantly 
overestimates the solid energy yield and densification ratio given 
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Fig. 1. Predicted solid energy yield (-) versus mass loss (-) in (a). Energy 
densification ratio (-) versus mass loss (-) in (b). Experimental data from (Almeida 
et al„ 2010; Arias et al„ 2008; Bergman and Kiel, 2005; Medic et al„ 2012; Ohliger 
et al„ 2012; Prins, 2005). 



Fig. 2. Solid and volatile energy yield (left axis) and heating value of volatiles 

the uncertainties associated with the experimental data and mod¬ 
eling approach. 

One reason the present model predictions may over-estimate 
the energy yield of the solid product is related to the fitted compo¬ 
sitions of VI and V2. As was noted in (Bates and Ghoniem, 2012) 
the fitted volatile composition predicts an evolution in solid 
product composition which slightly over-estimates the carbon 


yield in the solid product particularly at higher levels of mass loss. 
Since the heating value of the product is strongly correlated to its 
carbon content, over-predicting the carbon yield also over-predicts 
the energy yield. Additionally, differences between the model and 
experimental data can be explained by differences in the feedstock. 
The volatile composition data used in this model was fitted for data 
from willow whereas the experimental data used for comparison 
includes a variety of feedstocks such as beech, corn stover, euca¬ 
lyptus, and willow. 

At very high levels of mass loss (40-50%), which take 40 min 
under these kinetically controlled conditions, the model predic¬ 
tions appear to agree with the data within 1%. Overall, the agree¬ 
ment over the range of mass loss suggests the model can 
satisfactorily predict the evolution of the energy yield of the solid 
product during willow torrefaction. 

3.2. Volatile product energy yield 

The energy yield, r/vtot, (left axis) and average heating value, 
HHVvtot. (tight axis) of the cumulative volatile products as a func¬ 
tion of mass loss, VW are shown in Fig. 2 assuming the same ther¬ 
mal conditions as those in Fig. 1. Initially, the released volatiles 
consist only of VI (heating value of 4.43 MJ/kg). With higher levels 
of mass loss, increasing amounts of V2 (heating value of 16.3 MJ/ 
kg) are liberated, and thus the average heating value of the cumu¬ 
lative volatiles increases. For mass loss between 0 and 50%, the 
average heating value of the total volatiles ranges between 4.43 
and 10.6 MJ/kg. These model predictions for the average heating 
value of the volatiles compares very well with the experimental 
estimates made in (van der Stelt, 2011) and (Prins, 2005) where 
the LHV of volatiles produced during beech and willow torrefaction 
was estimated/measured to be between 1-8 and 4.9-10.6 MJ/kg, 
respectively (depending on the method of calculation and temper¬ 
ature conditions). (Neves et al„ 2011) reported a LHV between 2 
and 4 MJ/kg for volatiles released at a peak pyrolysis temperature 
of 300 °C. The present modeling results also agree with the exper¬ 
imentally observed trend that increased levels of mass loss are 
associated with increased volatile heating value. 

Fig. 2 also depicts the balance of energy in the solid and volatile 
products. As torrefaction proceeds with increasing mass loss, a 
higher fraction of the initial energy (i.e. heating value) of the solids 
is contained in the volatiles. After mass loss of 50%, the volatiles 
contain a significant percentage (28%) of that contained in the ini¬ 
tial solid product energy (18.7 MJ/kg). A mass loss fraction of less 
than 30% limits the fraction of energy contained in the volatiles 
to between less than 12% of the initial heating value of the solid 
product. 

3.3. Thermochemistry 

The results of thermochemistry for each reaction step (one and 
two) as a function of temperature are discussed in Section 3.3.1. 
Then, model results for the overall magnitude and rates of heat re¬ 
lease are discussed in Section 3.3.2. 

3.3.1. Enthalpy of reactions under isothermal conditions 

The predicted enthalpy of reaction for each step (based on Eqs. 
(18) and (19) and at the final char yield (based on Eq. (20)) are 
shown in Fig. 3a-c. 

AH ril calculated from Eq. (18) and plotted in Fig. 3a is the 
reaction enthalpy associated with the more rapid first stage 
(A -* B + V 1) corresponding to the decomposition of the hemicellu- 
lose fraction. Using the Boie correlation (HHV Eq. (1) from Table 3) 
the model predicts that depending on the temperature, this stage is 
exothermic (-40 to -180kJ/kg A ). This is in good quantitative 
agreement with the model predictions using the Friedl et al. corre- 
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Fig. 3. Enthalpy of reaction (inj kg ’) at temperature (°C) for torrefaction steps one 
(a), two (b), and at completion (c). 

lation (HHV Eq. (2) from Table 3) showing that the first stage is 
slightly exothermic (-50 to -92 kJ/kg A ). Applying the IGT correla¬ 
tion, (HHV Eq. (3) from Table 3) it is predicted that the first stage 
is exothermic (-65 to -280 kJ/kg A ). The exothermicity of the first 
stage can be explained by the composition of the volatile products 
released during this state. The first stage comprises of dehydration 
and decarboxylation reactions associated with the decomposition 
of the hemicellulose and (to a lesser extent) the cellulose fraction 
contained in the willow (Bates and Ghoniem, 2012; Chen and 
Kuo, 2011). As a result, the volatiles (VI) are comprised primarily 
of water, acetic acid, and carbon dioxide which have relatively 
low heat of formation (see Table Al). The predicted exothermicity 
of this first stage is consistent with the experimentally measured 
exothermicity of pure hemicellulose (xylan) pyrolysis bewteen 
200-300 °C (Khezami et al., 2005; Yang et al., 2007). 

AH r , 2 given by Eq. (19) and plotted in Fig. 3b is the reaction en¬ 
thalpy associated with the slower second stage (B -► C+ V2), corre¬ 
sponding to the decomposition of the remaining cellulose. Using 
the Boie correlation the model predicts that the second stage is sig¬ 
nificantly endothermic (180 to 206 l<J/kg B ) and has minimal tem¬ 
perature dependence. On the other hand, when applying the 
Friedl et al. correlation, the model predicts a shift from endother¬ 
mic to exothermic behavior at 245 °C. The predicted enthalpy of 
reaction ranges from 350kJ/kg B at 200 °C to -630 kJ/kg B at 
300 °C. Moreover, model predictions using the IGT correlation also 
show an exothermic second stage between -125 and -155 kJ/kg B . 
Thus, the application of various HHV correlations does not provide 
consensus on whether the second stage is either endothermic or 


exothermic. Pure cellulose undergoes competing and overlapping 
endothermic volatile formation and exothermic char formation 
reactions (Cho et al., 2010; Milosavljevic et al„ 1996; Mok and 
Antal, 1983). As a result, the overall enthalpy of cellulose pyrolysis 
also ranges from exothermic to endothermic (-500 to 800 kj/kg) 
depending linearly on the char yield (Cho et al., 2010). Therefore, 
while the model predictions for AH r , 2 are sensitive to the assumed 
HHV correlation they are reasonably within the range of experi¬ 
mentally observed data for cellulose pyrolysis thermochemistry. 

AH rfinal calculated from Eq. (20) and plotted in Fig. 3c) is the 
cumulative reaction enthalpy (per kg of initial biomass) occurring 
at the completion of both stages (A -► VI + C + V2) and correspond¬ 
ing to complete decomposition of the hemicelullose and cellulose 
fractions. Using the Boie correlation the model predicted that the 
enthalpy of torrefaction at completion ranges between 150 and 
-56 kj/kg for torrefaction temperatures of 200 and 300 °C, respec¬ 
tively. On the other hand, using the Friedl et al. correlation pre¬ 
dicted greater temperature dependence, varying between 275 kj/ 
kg endothermic at 200 °C to -540 kj/kg exothermic at 300 °C. 
Applying the IGT correlation, the model predicted the enthalpy of 
reaction at completion varying from -182 to -387 kj/kg. There is 
an almost identical variation with temperature when the Boie 
and IGT correlations are applied: in both cases the heat of reaction 
at completion becomes increasingly exothermic with temperature 
at a rate of -2.1 kJ/kg-K. Meanwhile, applying the Friedl et al. cor¬ 
relation predicts that temperature has a much stronger effect on 
the enthalpy of reaction (-8.2 kJ/kg-K). This might be explained 
by the fact that the Friedl correlation was fitted for raw biomass 
samples and is therefore perhaps not valid for estimating the heat¬ 
ing values of char-like materials. For example, with a char compo¬ 
sition of Yj,char = 0.92, 0.01, 0.016,0.01, 0.043 (/ = C, H, O, N, Ash) the 
Friedl et al. correlation estimates a heating value of 31.9 ± 0.33 MJ/ 
kg which is significantly lower than that estimated by the Boie cor¬ 
relation 33.3 ± 0.73 MJ/kg. 

Overall, when applying the three HHV correlations, the model 
predictions agree on the trend that both AH M and AH r , fina i, become 
increasingly exothermic with higher temperatures. The tempera¬ 
ture dependence is attributable to three factors: First, the two- 
stage competitive reaction mechanism results in a strongly tem¬ 
perature-dependent product distribution described by the instan¬ 
taneous fractional yields, /?, v, y, {, which are terms in Eqs. (18)- 
(20). As was discussed in (Bates and Ghoniem, 2012), higher tem¬ 
peratures shift the product distribution towards volatile products, 
a trend which causes increased exothermicity. Second, the compo¬ 
sition of the pseudo components B and C (which is accounted for in 
H b and H c ) in Eqs. (18)-(20) depend on /?, v, y, f. At higher temper¬ 
atures, B and C show increased carbon content and therefore a 
higher heat of formation. Lastly, the sensible enthalpy component, 
contained in the terms H A , H s , H c , H vu H v2 , is inherently tempera¬ 
ture dependent. This was shown by (Haseli et al., 2011) to signifi¬ 
cantly affect the heat of pyrolysis at elevated temperatures. Since 
the temperature range is relatively small (200-300 °C), this factor 
is less significant than the two aforementioned factors. 

A couple of works in the literature have discussed the temper¬ 
ature dependence of the heat of torrefaction. van der Stelt (2011) 
and Ohliger et al. (2012) found results which agree with the pres¬ 
ent model predictions: higher temperatures were associated with 
more exothermic reactions. Both authors hypothesized that at 
higher temperatures, exothermic secondary reactions occur be¬ 
tween the volatiles and the char. Ohliger et al. (2012) cites results 
from (Rath et al., 2003) which studied non-isothermal (a tempera¬ 
ture ramp at 10°C/min) beech and spruce pyrolysis for tempera¬ 
tures between 100 and 500 °C. However, it is not clear whether 
these secondary reactions are actually significant under the much 
milder temperatures associated with torrefaction (200-300 °C). 
The present modeling results offer an explanation for the temper- 
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ature-dependence of the heat of torrefaction based on the evolu¬ 
tion of primary volatile products alone. 


3.3.2. Dynamics of heat release during torrefaction 
3.3.2.1. Cumulative heat release. Fig. 4a demonstrates the cumula¬ 
tive (integrated) heat release q r (in units of J per kg of initial solid 
mass) for non-isothermal torrefaction conditions lasting 60 min 
with a heating rate of 10 °C/min, initial particle temperature of 
200 °C, and final particle temperature of 280 and 300 °C. Over this 
period, a mass loss of 36% and 56% of the original solid mass occurs 
for 280 and 300 °C, respectively. For both thermal conditions (280, 
300 °C), the cumulative heat release is net exothermic for the en¬ 
tire hour. Using the Boie HFIV correlation, the model predicted a 
maximum in cumulative heat release. This maximum corresponds 
to the end of stage one (which is the exothermic stage) and the 
beginning of stage two (which is the endothermic step). On the 
other hand, applying the Friedl et al. and IGT HFIV correlations, 
the model predicts that the cumulative heat release monotonically 
increases because both stages are exothermic at these tempera¬ 
tures. Depending on the assumed HHV correlation, by the end of 
the 60 min, the cumulative heat release at 300 °C is 65, 325, and 
418 kj/kg Boie and Friedl et al. and IGT correlations, respectively. 

The cumulative heat release q r is plotted versus the solid mass 
loss in Fig. 4b under the same thermal conditions as those in 
Fig. 4a. Depending on the final temperature, 280 or 300 °C, the 
reaction time required to achieve a mass loss of 40% is 77.5 and 


27.5 min under these kinetically controlled conditions, respec¬ 
tively. The model results demonstrate that for an equivalent level 
of mass loss, a higher torrefaction temperature results in a slightly 
higher cumulative heat release. Using the Boie HHV correlation, 
the model predicted a maximum exothermic point of torrefaction 
occurs for a mass loss between 27% and 32%. The Friedl et al. and 
IGT HHV correlations do not show maxima. Good quantitative 
agreement on the magnitude of cumulative heat release exists be¬ 
tween the predictions from the three HHV correlations over this 
mass loss range (0-40%). The predicted heat release is between 
0-112, 0-251, and 0-276 kj/kg when applying the Boie and Friedl, 
and IGT HHV correlations, respectively. 

The range of predictions provided by the three correlations 
agrees well with data from (Ohliger et al., 2012) which estimated 
a cumulative heat of release between 50 and 200 kj/kg for a mass 
loss of 33% during beech torrefaction. The model predictions also 
agree well with the exothermic cumulative heat of torrefaction 
estimate by (van der Stelt, 2011) from a single particle beech wood 
torrefaction of 0-200 kj/kg. The modeling results clearly demon¬ 
strate that increasing the torrefaction severity (as quantified by 
mass loss) results in greater cumulative heat release- lending sup¬ 
port to one of the main conclusions made in (Ohliger et al., 2012). 


3.3.22. Heat release rate. Fig. 5a-b illustrate the instantaneous heat 
release rate, (W/kg), versus mass loss and time, respectively, 
under the same non-isothermal conditions as those in Fig. 4. The 


I 



Time (seconds) 



Mass loss (Y| OSS ) 


Fig. 4. Cumulative heat release (J/kg initial solid) versus time time (a) and mass loss (b). T initia , = 200 °C, Heating rate = 10 °C/min. 7), = 280, 300 °C shown as green and red 
lines, respectively. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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heat release rate reaches a peak as the particle approaches the final 
temperature (at 8 and 10 min for 280 and 300 °C, respectively). 
Afterwards, the particle is maintained at its final temperature, 
and ^ decreases with time past this point. When the Boie HHV 
correlation is used, the model predicts that the heat release rate 
becomes negative as it enters the second endothermic stage. For 
a final torrefaction temperature of 280 °C, applying the Friedl 
et al., Boie, or IGT HHV correlations, the model predicts peak heat 
release rates of 170, 200, and 330 W/kg, respectively. At a higher 
final temperature of 300 °C, the predicted peak heat release rates 
are higher at 274, 327, and 565 W/kg, respectively. This clearly 
demonstrates that raising the peak torrefaction temperature in¬ 
creases the peak heat release rate. Analysis of Eq. (15) explains this 
trend: a higher final torrefaction temperature causes more rapid 
mass loss kinetics (i.e. reaction rate). Since the first stage is exo¬ 
thermic (AH r>1 < 0), a higher reaction rate results in a more rapid 
heat release rate. 

Since the peak heat release rate increases with particle temper¬ 
ature, the exothermic reactions predicted by the model could be 
expected to significantly affect the particle temperature especially 


in the presence of heat transfer limitations. Based on a specific heat 
capacity of raw wood (1500-2300 J kg 1 K '), a local exothermic 
heat release of 140 kj/kg would cause a temperature rise of up to 
60-90 Kelvin. In order to fully quantify the effect of this heat re¬ 
lease, a more detailed single particle model is needed to account 
for these transport limitations. A follow up publication will be de¬ 
voted to describing these coupled phenomena. 

3.4. Model evaluation and limitations 

The present model provides a method to predict the energy 
yield of the solid and volatile willow torrefaction products based 
on fitted volatile composition and mass loss kinetics parameters. 
It confirms the overall exothermicity of willow torrefaction and 
provides a fundamental explanation for the temperature-depen¬ 
dence of the reaction thermochemistry. Model predictions for the 
solid energy yield, rj s and energy densification ratio, R H hv compare 
satisfactorily with the experimental data even for feedstocks other 
than willow. Predictions for the volatile energy yield, rjy tot and 
cumulative heat release, q r , when using either of the three HHV 
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Fig. 5. Heat release rate (W/kg initial solid) versus mass loss (a) and time (b). Positive denotes exothermic. T iniHal = 200 °C, Heating rate = 10 “C/min. T, = 280,300 °C, shown in 
green and red lines, respectively. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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correlations appear to agree with the available literature data from 
beech and willow torrefaction. The predicted reaction enthalpies 
(AH rl , AH ri2 , AH rfinal ). are within the range of published data for 
torrefaction and pyrolysis of woody feedstocks summarized in 
Table 1. 

However, the exothermicity/endothermicity of the second reac¬ 
tion AH r>2 and- as a result- the rate of heat release, ^ appear sen¬ 
sitive to the assumed HHV correlation. This highlights the 
importance of two topics for future research. First, further experi¬ 
mental works applying differential thermal analysis (DTA) or dif¬ 
ferential scanning calorimetry (DSC) under torrefaction 
conditions to a variety of biomass feedstocks are needed to enable 
verification of these heat release rate predictions. Second, while 
correlations relating the ultimate analysis of a solid fuel to its heat¬ 
ing value have been developed for coal or raw biomass, none exist 
specifically for torrefied biomass. These HHV correlations are 
important because the errors in the HHV estimate are directly 
and linearly related to the errors in the estimate of the enthalpy 
of formation of the solid products. Since the magnitude of heat re¬ 
lease is less than 4% of the HHV of the dry willow (—18.5 MJ/kg), 
the development of improved correlations will enable more accu¬ 
rate thermochemical modeling of the torrefaction process. 

4. Conclusion 

The development of a thermochemical model coupled to an 
existing model for the evolution of volatile and solid matter during 
willow torrefaction was presented. Composition dependent ther¬ 
modynamic properties of the torrefaction pseudo components 
were estimated, enabling a description of the rates/magnitude of 
heat release. The model enables a more robust energy accounting 
for multi-scale torrefaction models (e.g. single particle, reactor, 
flowsheet models). 

The model predicts an exothermic heat release in the first stage, 
between 40 and 280 kJ/kg biomass . The exothermicity of the second 
stage depends on temperature and the assumed heating value cor¬ 
relation. Peak heat release rate increases with higher torrefaction 
temperatures. 
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